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The energy spectrum for three species of identical fermionic atoms close to a Feshbach resonance 
is computed by use of a nonperturbative flow equation. Already a simple truncation shows that for 
large scattering length \a\ the lowest energy state is a "trion" (or trimer) bound state of three atoms. 
At the location of the resonance, for \a\ oo, we find an infinite set of trimer bound states, with 
exponentially decreasing binding energy. This feature was pointed out by Efimov. It arises from 
limit cycle scaling, which also leads to a periodic dependence of the three body scattering coupling 
on In \a\. Extending our findings by continuity to nonzero density and temperature we find that a 
"trion phase" separates a BEC and a BCS phase, with interesting quantum phase transitions for 
r = 0. 



I. INTRODUCTION 

Ultracold fermion gases with three hyperfine species 
near a Feshbach resonance are far from being understood 
as well as the analogous Bareen-Cooper-SchriefFer (BCS) 
- Bose-Einstein condensate (BEC) crossover with only 
two species [ij. Only recently such systems were investi- 
gated experimentally, showing interesting phenomena be- 
yond the two species [21 . It is a challenge for both theory 
and experiment to understand the phase diagram of such 
systems. Although the physics of three-body scattering 
can be described with different methods reaching from 
quantum mechanics to quantum field theory, it is difh- 
cult to account for many-body effects in thermodynamic 
systems with a nonzero density n and temperature T. 
For ultracold fermions in an optical lattice a phase with 
"trions" as well as a phase with "color superfluidity" was 
found in a variational approach Q , see also Q • 

The functional renormalization group provides a pow- 
erful method with the potential to describe both the scat- 
tering amplitudes and binding energies of the few-body 
sector and the many-body physics at nonzero density and 
temperature. A nonperturbative exact flow equation in- 
cludes the effect of quantum and statistical fluctuations 
by following the evolution of the average action Tk from 
"microphysics" to observable "macrophysics" [Sj]. 

In this paper we concentrate on the functional integral 
description of energy levels of the dimer and trimer bound 
states ( "vacuum physics" ) . We demonstrate that a rather 
simple truncation of the flow equation can describe the 
dependence of various energies on the scattering length 
a, as shown in Fig. [T] At the "unitarity limit" a — > cxd 
for B = Bq we also find an infinite tower of trimer bound 
states, with energies given by 



En = exp ( — — ) Eq, 



(1) 



done within the same framework of truncations, simply 
by changing the chemical potential and the temperature. 

Already with the present "vacuum computation" we 
can conclude by continuity that for small temperature 
T two phases with "extended superfluidity" are sepa- 
rated by a "trion phase" with a different symmetry. The 
extended superfluidity is characterized by several Gold- 
stone bosons, corresponding to the symmetry breaking 
SU(3) X U(l) SU(2) X U(l) and is analogous to the 
BEC and BCS superfluidity [1]. 
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the spectrum found by Efimov in a quantum mechani- 
cal computation @, for a review see Q. The functional 
renormalization group is thus a promising tool for an 
investigation of the many-body physics, which can be 
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FIG. 1: (Color online) Energy spectrum for three fermion 
species close to a common Feshbach resonance. The solid 
line gives the energy per atom of the trion with respect to 
the fundamental fermion (zero energy line). The dashed line 
is the energy of the dimer per atom. Here ao is the Bohr 
radius which we also use as the ultraviolet scale A = I/uq. 
Parameters are chosen to correspond to ^Li in the {mp ~ 1/2, 
mp = — l/2)-channel [^. 



II. METHOD AND APPROXIMATION SCHEME 

In a formulation with scale dependent fields, the effec- 
tive average action (also called flowing action or running 
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action) obeys the flow equation 

Here, the (generalized) spinor $fc describes afl fields in 
the problem and STr sums over momentum q and Mat- 
subara frequency (jo as well as the internal degrees of 
freedom such as species of fields, with a minus sign for 
fermions. By Rk we denote an infrared regulator which 
suppresses fluctuations with momenta k^- At the 

ultraviolet scale fc = A the average action equals the mi- 
croscopic action Fa — S. In our case this momentum 
scale is of the order of the inverse Bohr radius A — I/oq, 
ao ~ 5 X lO^^^m. Including more and more fluctuation 
effects by flowing to smaller cutoff scales k, we finally 
arrive at the full quantum effective action T = Tk=o, 
which is the generating functional for the one-particle ir- 
reducible correlation functions. Precise knowledge of T 
corresponds to a "solution" of the theory. Finding ex- 
act solutions to the flow equation ^ is not possible in 
practice but one can find approximate solutions by trun- 
cating the space of functionals F^ to a manageable size. 
Reviews on the method of functional renormalization can 
be found in 

We propose a scheme to solve Eq. ([2]) approximately 
for the case of three species of identical fermions. Apply- 
ing it to the vacuum, where n = T = 0, yields a simple 
and natural description of the scattering properties. We 
use the following truncation for the average action 

+he,,k (0*^jVfe - M*rk) /2 + (0t0)' /2 
+X* {dr - A/3 + m\) x + g {(t)li^*iX - '/'^V'^X*) 

+A^^(0*^,r(/.,^,)}- (3) 

Here we use natural nonrelativistic units with h ^ ks ^ 
2M = 1, where M is the mass of the original fermions. 
The integral in Eq. ([3]) goes over homogeneous space 
and over imaginary time as appropriate for the Matsub- 
ara formalism = / (fix Jq^'^ dr. The symbol A is the 
Laplace operator. We denote the (Grassmann) field de- 
scribing the three fermion species hy ip = (■01, '02, "03): 
the field of bosonic bound states by = (0i,02,03) = 
('02'03j V'3'01, V'iV'2) and also include a single component 
fermionic field x- This "trion" field represents the totally 
antisymmetric combination '0i02'03- 

On the level of the three-body sector, the symmetry 
of the problem would allow also for a term ~ ■0'''00t0 in 
Eq. ([3]). This term plays a similar role as for the case 
of two fermion species, where it was investigated in p^ . 
We expect that the qualitative features of the three-body 
scattering are dominated by the term A^^ in Eq. ([3]). 



The quantitative influence of a term ^ ip^ ilj(f>^ (p on the 
flow equations will be investigated in future work. 

We assume in this paper that the fermions 0i , 02 , and 
■03 have equal mass M and chemical potential /i. We 
also assume that the interactions are independent of the 
spin (or hyperspin) so that our microscopic model is in- 
variant under a global SU(3) symmetry transforming the 
fermion species into each other. While the fermion field 
■0 = (■01,^02, ^03) transforms as a triplet 3, the boson field 
— (01, 02, 03) transforms as a conjugate triplet 3. The 
trion field x is a singlet under SU(3). In concrete experi- 
ments, for example with ^Li [2], the SU(3) symmetry may 
be broken explicitly since the Feshbach resonances of the 
different channels occur for different magnetic field values 
and have different widths. For nonzero density this pat- 
tern may lead to interesting phase diagrams with various 
types of superfluidity . 

In addition to the SU(3) spin symmetry our model is 
also invariant under a global U(l) symmetry ip —^ e*"^0, 
e^*"0, and x ~^ e'^*"x. The conserved charge related 
to this symmetry is the total particle number. Since we 
do not expect any anomalies the quantum effective action 
F = Ffe=o will also be invariant under SU(3) x U(l). 

Apart from the terms quadratic in the fields that de- 
termine the propagators, Eq. ^ contains an interaction 
between bosons ~ A^, a term that describes scattering 
between bosons and fermions Xff,^ and the Yukawa- 
type interactions ~ h and ^ g. The energy gap parame- 
ters for the bosons and for the trions are some- 
times written as = Vcj, — 2/i, = i/-^ — 3/i, showing 
an explicit dependence on the chemical potential /i. 

In Eq. ([3]), the fermion field x can be "integrated 
out" by inserting the (0, 0)-dependent solution of its 
field equation into F^:. For — 00 this results in a 
contribution to a local three-body interaction, shifting 
A^^ Xcj,,p — g'^/m'^. Furthermore for A^ = one may 
integrate out the boson field 0, such that (for large mp 
one replaces the parts containing and x i^i by an 
effective pointlike fermionic interaction 

rfc,int=^^A^(V'V)' + ^A3(V'V)', (4) 

with 




We note that the contribution of trion exchange to A^^ 
or A3 depends only on the combination g^/m^. The sign 
of g can be changed by x ^ — x, and the sign of g^ can 
be reversed by a sign flip of the term quadratic in x- 
Keeping the possible reinterpretation by this mapping in 
mind, we will formally also admit negative g^ (imaginary 

g)- 

At the microscopic scale A; = A, we use the initial val- 
ues of the couplings in Eq. © 5 = A^ = A^^ = and 
TO^ —> 00. Then the fermionic field x decouples from the 
other fields and is only an auxiliary field which is not 
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propagating. However, depending on the parameters of 
our model we will find that x, which describes a compos- 
ite bound state of three original fermions x = V'iV'2'03j 
becomes a propagating degree of freedom in the infrared. 
The field x is called trion. The initial values of the boson 
energy gap and the Yukawa coupling h will determine 
the scattering length a between fermions and the width 
of the resonance, see below. The pointlike limit (broad 
resonance) corresponds to i^^ ^ oo, /i^ — s- oo where 
the limits are taken such that the effective renormalized 
four fermion interaction remains fixed. In Eq. ([3]) we 

use renormalized fields 4> — ^^'^{k)<j}-, "0 = ^J,'^' 



X = A'^'\k)x, with A^{A) = A^{A) = 4(A) = 1, and 



renormalized couplings mi — mi/A^, h = h/{A^/'^A^), 



A^^ = X^^/{Arj,A^) 



rnl/A^, g 



-gl^A^rATA^r)^ and 



To derive the flow equations for the couplings in Eq. 
^ we have to specify an infrared regulator function i?^ . 
Here we use the particularly simple function i?^ = r(fc^ — 
p^)0(fc^— p^), where r = 1 for the fermions ip, r = 1/2 for 
the bosons </>, and r = 1/3 for the composite fermionic 
field X- This choice has the advantage that we can derive 
analytic expressions for the flow equations and that it 
is optimized in the sense of We use the method 

of "(re-)bosonization" neglecting correction terms 

which have only a minor quantitative effect. 



III. FLOW EQUATIONS FOR TWO-BODY 
SECTOR 



We can now insert our ansatz equation ^ into the flow 
equation Functional derivatives with respect to the 
fields for zero temperature T — and density n = lead 
us to a system of ordinary coupled nonlinear differential 
equations for the couplings m^, h, A^, A^, A^, fn^, g 
and A^^. One finds that the propagator of the original 
fermions ip is not renormalized, A^{k) = 1, h{k) = h. 
The flow equations for the couplings determining the two 
body sector, namely the boson gap parameter (with t = 
ln(fc/A)) 



decouple from the other flow equations, 
solved analytically 

' 67r2 



They can be 




Using this explicit solution ^ we can relate the initial 
value of the boson gap parameter as well as the 
Yukawa coupling h to physical observables. The inter- 
action between fermions is mediated by the exchange 
of the bound state (f). Again in terms of bare quanti- 
ties, the scattering length between fermions is given by 
a — —h?/ (Sttto^) where the coupHngs h and are eval- 
uated at the macroscopic scale fc = and for vanishing 
chemical potential /i = . This fixes the initial value 



m^(A) = -^a-i 



—a + —A - 2^. 



(9) 



Stt ' 67r2 

In addition to the /i-independent part we have added here 
the chemical potential term — 2/i where the factor 2 ac- 
counts for the bosons consisting of two fermions. In vac- 
uum, the gap parameter of the bosons is proportional 
to the detuning of the magnetic field m^(fc = 0, ^ = 0) = 
Plm{B — Bq). Here hm is the magnetic moment of the 
bosonic dimer and Bq is the magnetic field at the reso- 
nance. From 

1 

"""8^MA/(S-So) 

one can read off that h? is proportional to the width of 
the resonance. We have now fixed all initial values of 
the couplings at the scale fc = A or, in other words, the 
parameters of our microscopic model. 

At vanishing density n = 0, the chemical potential is 
negative or zero, < 0, and will be adjusted such that 
the lowest excitation of the vacuum is a gapless propa- 
gating degree of freedom in the infrared, i. e. at fc = 0. 
Depending on the value of this lowest energy level 
may be the original fermion -0, the boson 0, or the com- 
posite fermion x- 



_ ^ (f!\ IV. THREE-BODY PROBLEM 

- 6^(fc2-^)2' 

Now that we have solved the equations for the two- 
body problem within our approximation, we can address 
the three-body sector. It is described by the flow equa- 
tions for the trion gap parameter, 

^7 P fc5 p. 2 6g^ fc^ 2 ni^ 

^*^^ = -6^(fc^ ^*"^x = ^(3P_2^_,2..2^)2+%-x^ (11) 



and the boson wave function renormalization 



4 



and the Yukawa-type coupling, 

dtg = (770 + 7]x)| + "^^9ta 



6fc2 - 5/^ + 2ml 



(12) 



These equations are supplemented by the anomalous di- 
mension 

24^2 



A. 



X 



and the variable 



7r2 (3/c2 - 2/1 + 2m2 )3 
/fc5 9fc2 - 7/i + 4m2 



127r25 (fc2 - /i)3 (3/c2 _ 2/t + 2m2 )2 ' 



(13) 



(14) 



Here a determines the scale dependence of the trion field 

X 



dtX = {dta)(l),^pi + 77xX/2. 



(15) 

While the second term in Eq. psp is the usual wave- 
function renormalization, the first term describes a non- 
linear change of variables. It is chosen such that the flow 
of X^jj, vanishes on all scales, Xcj,^{k) = 0. The scatter- 
ing between bosons and fermions is then described by the 
exchange of the trion bound state x- This reparametriza- 
tion, which is analogous to rebosonization (ill . [l7j , is cru- 
cial for our description of the system in terms of the com- 
posite trion field x- Note that the flow equations pT|) . 
P^ . p^ . that describe the three-body sector are 
independent from the flow of the boson-boson interaction 
that belongs to the four-body sector. 
Since the three-body sector is driven by fermionic and 
bosonic fluctuations, it is not possible to find simple an- 
alytic solutions to the flow equations in the general case. 
However, it is no problem to solve them numerically. This 
is most conveniently done using again "bare" couplings. 
As a general feature we note that a negative chemical 
potential /i acts as an infrared cutoff for the fermionic 
fluctuations while a positive value of the bosonic gap m'^ 
suppresses bosonic fluctuations in the infrared. We can 
find numerical solutions for different scattering length be- 
tween the fermions a and varying chemical potential /i. 
We use an iteration process to determine the chemical po- 
tential fj- < where the lowest excitation of the vacuum 
is gapless. On the far BCS side for small and negative 
scattering length a where this lowest excitation 

is the fermion ip^ this implies simply /i = 0. On the far 
BEG side for small and positive scattering length a — > 0+ 
the lowest excitation is the boson cf), implying — 0. 
We can then use our analytic solution of the two-body 
sector Eq. ([5]) to obtain the chemical potential /i from 
the condition m| = at the macroscpic scale fc = 0. In 
the limit A/|/i| — > 00 we find 



/l4 



l/l2 



(32^)2 



167r 



(16) 



In the broad resonance limit ft.^ 
the well-known result /i = — 1 /a2 . 



00 this reduces to 
For scattering lengths 



close to the Feshbach resonance a,. 



< a < 



find that the lowest vacuum excitation is the trion x- 
Numerically, we determine fi from the implicit equation 
m2 = at fc = 0. Our result for /i obtained for = 100 A 
as a function of the inverse scattering length is shown 
in Fig. [2] With the choice A — 1/ao, this value of 
corresponds to the width of the Feshbach resonance of 
^Li atoms in the (m^ = 1/2, mp = — l/2)-channel 
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FIG. 2: (Color online) Dimensionless chemical potential 
— y'JJI\/A in vacuum as a function of the dimensionless scat- 
tering length a^^/A. For comparison, we also plot the result 
for two fermion species where the original fermions are the 
propagating particles for < (dotted) and composite 
bosons for > (dashed). The inset is a magnification of 
the little box and shows the energy of the first excited Efimov 
state. 



For small Yukawa couplings h^/A^ 1, or narrow res- 
onances we find that the range of scattering length where 
the trimer is the lowest excitation of the vacuum in- 
creases linear with Ji^ [l^, |T9]. More explicit, we find 



'■cl 



-0.0015 /i2, a 



c2 



0.0079/1 . However, for very 
broad Feshbach resonances /A^ 1 the range depends 
on the ultraviolet scale a^/, ~ ^- We show this be- 
havior in Fig. [3] 

The chemical potential at the unitarity point a^^ — 
is plotted in the inset of Fig. O It increases with the 
width of the resonance similar to /iy = —3.5 x lO^^ft.* for 
small h? I A. It also approaches a constant value which 
depends on the cutoff scale ^jj ^ A'^ in the broad reso- 
nance limit 00. 

The dispersion relation for the atoms, dimers and 
trimers can be computed by analytical continuation, 
T = it, of the inverse propagators to "real frequencies" 
Lu. In our truncation, they read (in terms of renormalized 
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3: (Color online) Interval of scattering length < 
< a~2 where the lowest vacuum excitation is the trimer 
fermion x (solid lines). For comparison we also plot the linear 
fits a~i^ = -0.0015/1^ and a'^^ = 0.0079 /i^ (dotted). The 
inset shows the chemical potential fiu/^'^ at the resonance 
= as a function of ft^/A in the same range of Ji^ /A. Here 
we also plot the curve /xu = —3.5 x 10~®/i* for comparison 
(dotted). 



fields) 

= 



Px = 



2M 
AM 
6M 



- /I, 



(17) 



We note that the functions i^^ and depend only on 
the particular combinations of /j,, lo, and jfi given above. 
This is a result of the symmetries of our problem. The 
real-time microscopic action S — Vk=h with t ~ —it 
is invariant under the time-dependent U(l) symmetry 



transformation of the fields ip 



iEt 



A2Et, 



and X ~^ e x if we also change the chemical poten- 
tial according to ^—ffl = fj, + E. Since the microscopic 
action has this symmetry and since we do not expect 
any anomalies, the quantum effective action T — Tk=o 
is also invariant under these transformations. In conse- 
quence, any dependence on /i must be accompanied by 
a corresponding frequency dependence, such that only 
the invariant combinations uj^ + ^J■, i^<t> + 2/i, lu^ + 3/i 
appear in the effective action. The relation between 
the dependence on u and on jr^ is a consequence of 
Galilean symmetry. Even though we compute directly 
with the flow equations only j/0(/i,u;0 = 0,p*^ = 0) and 
v^{^,ijj^ = 0,p^ = 0), we can use the symmetry informa- 
tion ( "Ward identities" ) for an extrapolation to arbitrary 
uj and . The dispersion relation for the bosons ^^(p^) 
follows from Eq. p?|) by solving P^{uj^ = E^) = and 
similarly the dispersion relation of the fermions £'^(p*^) 
and the trions £^{p^). 

We are interested in the dimer and trimer energy levels 



relative to the energy of the atoms. For this purpose we 
consider their dispersion relations at rest {ff = 0) and 
substract twice or three times the atom energy £^ = — ^, 



E^ = S^if = 0) + 3^i = ,y^{E^/3). 

When the dimers are the lowest states one has 
and therefore — 2fi — E^ 
~ implies ~ = E- 



(18) 

while for lowest trimers 
We have shown E^ and 
E^ in Fig. [TJ A typical value for A is of the order of 
the inverse Bohr radius a^^ . We also mention that the 
energy levels have been computed in the absence of a mi- 
croscopic three-body interaction, A3 (A) = 0. The same 
computation may be repeated with nonzero A3 (A), and 
thus microscopic parameters may be fixed by a compar- 
ison to the experimentally measured spectrum. 

So far we were only concerned with the lowest energy 
excitation of the vacuum. We found that near a Fesh- 
bach resonance this lowest excitation is given by a trimer. 
However, it is known from the work of Efimov I6l| that one 
can expect not only one trimer state close to resonance 
but a whole spectrum. In fact, the solution = 3fi 
may not be the only solution for the equation fixing the 
trimere energy levels, i. e. 



£x = ''x{f^ + l^x 



3/i. 



(19) 



For an investigation of the energy dependence of we 
use the symmetry transformation and shift the chemical 
potential hy fi = 11 + £^ = 0. We can therefore follow 
the flow at vanishing p, and Eq. ()19p turns into a simple 
implicite equation for the trion energy levels 



Ex = -x{^ 



E. 



(20) 



For fl — the flow equations simplify considerably. For 
example, Eq. ([8]) becomes 

(21) 



Mk) = 1 + ^ ( T - T 



We observe that the two-body sector is evaluated here for 
LUfj, = Utf, + 2^ = 0. The physical chemical potential re- 
mains, of course, negative in the vicinity of the Feshbach 
resonance, and our evaluation therefore corresponds to a 
positive energy Cij0 as compared to the lowest energy level 
of the trimer for which lu^ ^ 0. 

An especially interesting point in the spectrum is the 
unitarity limit, B — Bq, where the scattering length di- 
verges, = 0. At that point all length scales drop out 
of the problem and we expect a sort of scaling solution 
for the flow equations. In the limit /c — )■ the solution 
for A^{k) is dominated completely by the term with 1/k, 
A^{k) — h? / {Qir'^k), and we find vn?, — rf 
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For Ey^ the flow equations for the three-body sector 
simphfy and we find 



dt-ml 



dtg 



64 



_ H _ 



(22) 



In contrast, for ^ 0, as needed for Eq. (|^ . we 
will have to solve the flow with a nonvanishing "effec- 
tive chemical potential" fi — E-^/3, which will cause a 
departure from the scaling flow. 



smaller than For an evaluation of Eq. ([20)1 we may 

therefore use Eq. (p4|l with a specific value for fc, namely 
k'^ = —E^/3. The possible energy levels therefore obey 



Sk"^ + A-\k) mlik) = 



With 



dkA^ik) = 



24 fjk) 
1257r2 k^ 



(26) 



(27) 



one finds, up to oscillatory behavior an increase of 
fc-M . We can write Eq. in the form 



V. LIMIT CYCLE SCALING 

First, let us consider the scaling solution obeying Eq. 
. It is convenient to rescale the variables according to 
TO^ — fh'^{k/h)^, g"^ — g'^{k/hy^^ which gives the hnear 
differential equation 



5m ,2 



13 9 

' 25 ' ^ 



36 
25 



(23) 



The matrix on the right hand side of Eq. has the 

eigenvalues P1/2 = 9- (7/25)±«V4T9/25. The flow equa- 
tion leads therefore to an oscillating behavior. It is 
straightforward to solve Eq. explicitly. Restricting 
to real solutions and using initially (A) = wc find the 
following for the trimer gap parameter and coupling 



rhlik) 



f{k) = 



cos So In 




A 



(A) sin ( sq In 



(24) 



As it should be, the initial value to^(A) drops out of the 
ratio and we find for the three-body coupling 



A3(fc) 



4687r4 



sin (sq In 4 



(25) 

We obtain for the "frequency" sq = %/4T9/25 w 0.82. 
Since we use a truncation in the space of functionals 
this result is only a rough estimate. It has to be compared 
with the result of other methods which find sq ~ 1.00624 
[1) S HO] ■ Considered the simplicity of our approxima- 
tion, the agreement is quite reasonable. 

For a determination of the trion energy levels we have 
to solve the flow with an effective negative chemical po- 
tential jl = E^/'i. This acts as an infrared cutoff, such 
that the flow deviates from the limit cycle once fc^ « 
Qualitatively, the flow eventually stops once k becomes 



F{k) - 



cos I Sq In — 




(28) 



For small fc <C A we infer that F{k) has to vanish ~ 
— k/h?. Since F{k) is periodic, solutions will occur for 
roughly equidistant values in In ^ . 

For k <^h? the possible solutions simply correspond 
to F{k) = 0. The flrst solution with the largest k cor- 
responds to the ground state level with iJo < 0. The 
subsequent solutions obey 



So 



In- 



A 



In- 



A 



or 



E, 



n+l 



En 



exp 



27r\ 
— , En 

So , 



exp 



(29) 



27rn\ 

]Eo. 30 

So 



This corresponds to the tower of trimer bound states at 
the unitarity limit, with E^ approaching zero exponen- 
tially for n — > 00. For Eq. (|30|) we have actually taken 
into account all zeros of F{k). Since (f'{k) oscillates peri- 
odically, only half of these zeros correspond to g'^{k) > 0, 
while the other half has formally g^{k) < 0. We may 
use the mapping discussed after Eq. ^ to obtain an 
equivalent picture with positive g^. 

We may understand the repetition of states by the fol- 
lowing qualitative picture. The coupling to^, that is pro- 
portional to the energy gap of the trimer, starts on the 
ultraviolet scale with some positive value. The precise 
initial value is not important. The Yukawa- type coupling 
g vanishes initially so that the trion fleld x is simply an 
auxiliary fleld which decouples from the other flelds and 
is not propagating. However, quantum fluctuations lead 
to the emergence of a scattering amplitude between the 
original fermions ^ and the bosons (p. We describe this by 
the exchange of a composite fcrmion x- This leads to an 
increase of the coupling g"^ and a decrease of the trion gap 
rri^. At some scale ti — ln(fci/A) with kf « — /i the cou- 
pling crosses zero which indicates that a trion state x 
becomes the lowest energy excitation of the vacuum. In- 
deed, would we consider the flow without modifying the 
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m2(0/m2(0), r(i)/5'(0) 
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FIG. 4: (Color online) Limit cycle in the renormalization 
group flow at the unitarity point a^^ — 0, and for energy 
at the fermion threshold jl = fj, + E — 0. We plot the rescaled 
gap parameter of the trimer m^{t) (solid) and the rescaled 
Yukawa coupling g^{t) (dashed). The dotted curves would be 
obtained from naive continuation of the flow after the point 
where = 0. 

chemical potential, this would set an infrared cutoff that 
stops the flow at the scale ki ^ \/JJA ^''^'^ the trion x 
would be the gapless propagating particle while the orig- 
inal fermions ip and the bosons (j) are gapped since they 
have higher energy. 

Following the flow further to the infrared, we find that 
the Yukawa coupling decreases again until it reaches 
the point — at t — t[ (see Fig. HJ. Naive continua- 
tion of the flow below that scale would lead to < and 
therefore imaginary Yukawa coupling g. However, since 
the trion field x decouples from the other fields ioi g — 0, 
we are not forced to use the same field x before. We 
can simply use another auxiliary field X2 with very large 
gap = '^X2/^X2 to describe the scattering between 
fermions and bosons on scales t < t[. We are then in 
the same position as on the scale fc = A and the process 
repeats. Starting from a positive value, the rescaled gap 
parameter 771^2 decreases as the infrared cutoff k is low- 
ered. At the scale fc2 it crosses zero which indicates that 
there is a second trimer bound state in the spectrum with 
energy per original fermion E2 = — /i — fc|. Would we use 
the modified chemical potential ji = ^ + E2 = — fc|, the 
fiow would be stopped at the scale fc2 and the second 
trimer X2 would be the propagating degree of freedom. 
This cycle repeats and corresponds precisely to the limit 
cycle scahng of Eq. (^5]) . 

At the unitarity point with a^^ = and at the thresh- 
old energy E ^ —fi with // = /i + = 0, this limit cycle 
scaling is not stopped and leads to an infinite tower of 
trimer bound states. The energy of this states comes 
closer and closer to the fermion-boson threshold energy 
E = —fi. In our language we recover the effect first pre- 
dicted by Efimov p]. A similar limit cycle description 
of the Efimov effect for identical bosons was given in the 
context of effective field theory in [l^l ■ 



The infinite limit cycle scaling occurs only directly 
at the Feshbach resonance with = 0. Similar to a 
nonzero (modified) chemical potential jl, also a nonzero 
inverse scattering length provides an infrared cutoff 
that stops the fiow. For example for negative scattering 
length a and energy E with p, — /.i + E = 0, the solution 
of the two body sector Eq. ([5]) implies for the boson gap 
in the infrared = —Sira^^k/i + fc^ so that bosonic 
fiuctuations are suppressed in comparison to the unitar- 
ity point with = 0. This leads to a stop of the limit 
cycle scaling at the scale fc « 37r|a~^|/4 so that the spec- 
trum consists only of a finite number of trimer states. In 
the inset of Fig. [2] we show the numerical result for the 
modified chemical potential p, of the first excited Efimov 
trimer. 

VI. CONCLUSION 

Our simple truncation of an exact functional renor- 
malization group equation yields a quite detailed picture 
for the spectrum of bound states and the scattering of 
three species of identical fermions. Depending on the 
scattering length a, the lowest energy state may be a fun- 
damental fermion (atom) or a bosonic molecule (dimer) 
or even a composite fermion formed from three atoms 
(trimer). Trimers are stable (lowest energy state) for 

an interval {a^ci)~^ < < (ac2^)~^ which depends 
on the Yukawa or Feshbach coupling h related to the 
width of the Feshbach resonance. We have computed the 
binding energy Ei < of the trimer, cf. Fig. [1] For 

i'^'ci)~^ < < (0^2')"^ ^ second trimer bound state is 
present, with binding energy i?2. Close to the resonance, 
we find the (approximate) relations 

En+i = exp En, 

V So/ 

{a^Z'Y'-{ai%)-^e.p[-^y (31) 

Our truncation yields Sq — 0.82. This structure repeats 
for higher n, such that at the location of the Feshbach 
resonance for — > an infinite tower of trimer bound 
states builds up. The relations, Eq. (|3T|l. become exact 
for large n, if the quantity sq is computed precisely. 

The origin of this repetition of self-similar features on 
different scales is due to a limit cycle scaling behavior 
of the renormalization fiow. The same features appear 
periodically for sq In — -—nn, with fc the appropriate 
scale of the model. In other words, at the scale fcn+i = 

fc„ exp ^) one finds the same physics as at the scale 

fc„. This concerns all dimensionful quantities according 

to their dimension, e.g. En ~ fc^, ia''ci]2)~^ ~ ^n- ^Iso 
the three-body scattering length A3 depends periodically 
on the logarithm of the characteristic scale fc, with an 
overall factor fc""', cf. Eq. ([25]) . 

We have only considered the vacuum in the present 
paper, but we can already infer some interesting features 
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of what happens at nonzero density for the many-body 
ground state at T = 0. The quahtative properties for 
smah density follow from our computations by simple 
arguments of continuity. For small density, we find that 
a "trion phase" separates the BEC phase and the BCS 
phase. 

For increasing density, the chemical potential fi in- 
creases. The BEC phase occurs for small density for 
> (ac2)~^. Due to the symmetries of the micro- 
scopic action, the effective potential for the bosonic field 
depends only on the SU(3) x U(l) invariant combination 
p = (f)^(j) = -I- 02'/'2 + 03</>3- It reads for small /i — /io, 
with /xq the vacuum chemical potential, 



(32) 



where we use the fact that the term linear in p vanishes 
for ^ = poi i-e. ^^(/io) = 0. The minimum of the 
potential shows a nonzero condensate 



Po = 



■> U{po) = - — (m- Mo) , 



with density 



dUjpo) 
dp 



A, 



-{P - Mo)- 



(33) 



(34) 



For the BCS phase for a < Oci (od < 0), one has 
po = 0. Nonzero density corresponds to positive p — po, 
and in this region the renormalization flow drives al- 
ways to zero at some finite kc, with BCS spontaneous 
symmetry breaking (po > 0) induced by the flow for 
k < kc- Both the BEC and BCS phases are therefore 
characterized by superfluidity with a nonzero expecta- 
tion value of the boson field 0o 7^ with po — (f'o't'o- 

As an additional feature to the BCS-BEC crossover for 
a Fermi gas with two components, the expectation value 
for the bosonic field (po in the three component case also 
breaks the spin symmtry of the fermions (SU(3)). Due 
to the analogous in QCD this was called "color superflu- 
idity" [3. For any particular direction of (po a continu- 
ous symmetry SU(2) x U(l) remains. Acoording to the 
symmetry breaking SU(3) x U(l) ^ SU(2) x U(l), the 
effective potential has five flat directions. 

For two identical fermions the BEC and BCS phases 
are not separated, since in the vacuum either ^0 = or 
= 0. There is no phase transition, but rather a con- 
tinuous crossover. For three identical fermions, however, 

we find a new trion phase for (0^1^)^^ < a^^ < {o,^c2)~^- 
In this region the vacuum has po < and > 0. The 
atom fluctuations are cut off by the negative chemical po- 
tential and do not drive m| to zero, such that for small 
density remains positive. Adding a term m'^p to the 
effective potential (1321) we see that the minimum remains 



at Po = as long as > 2{p — po)- No condensate of 
bosons occurs. The BEC and BCS phases that show both 
extended superfluidity through a spontaneous breaking 
of the SU(3) X U(l) symmetry, are now separated by a 
phase where (po = 0, such that the SU(3) x U(l) symme- 
try remains unbroken (or will be only partially broken). 

Deep in the trion phase, e.g. for very small |a^^|, the 
atoms and dimers can be neglected at low density since 
they both have a gap. The thermodynamics at low den- 
sity and temperature is determined by a single species 
of fermions, the trions. In our approximation it is sim- 
ply given by a noninteracting Fermi gas, with fermion 
mass 3M and chemical potenial 3(/i — po)- Beyond our 
approximation, we expect that trion interactions are in- 
duced by the fluctuations. While local trion interactions 
~ (x*x)^ forbidden by Fermi statistics, momentum 
dependent interactions are allowed. These may, however, 
be "irrelevant interactions" at low density, since also the 
relevant momenta are small such that momentum depen- 
dent interactions will be suppressed. Even if attractive 
interactions would induce a di-trion condensate, this has 
atom number six and would therefore leave a Zg subgroup 
of the U(l) transformations unbroken, in contrast to the 
BEC and BCS phase where only Z2 remains. Further- 
more, the trions are SU(3)-singlets such that the SU(3) 
symmetry remains unbroken in the trion phase. The dif- 
ferent symmetry properties between the possible conden- 
sates guarantee true quantum phase transitions in the 
vicinity of Od and ac2 for small density and T = 0. We 
expect that this phase transition also extends to small 
nonzero temperature. 

While deep in the trion phase the only relevant scales 
are given by the density and temperature, and possi- 
bly the trion interaction, the situation becomes more 
complex close to quantum phase transition points. For 
a « flci we have to deal with a system of trions and 
atoms, while for a « ac2 a system of trions and dimers 
becomes relevant. The physics of these phase transitions 
may be complex and rather interesting. 

We finally comment on the precision of our compu- 
tation. At nonzero density, our truncation may be im- 
proved by including in the two-body sector a four fermion 
coupling ~ (V'^V')^: which may be partially bosonized in 
favor of a running h [2l| . Also the momentum depen- 
dence of the interactions ~ X^^ or A^ may be resolved 
beyond the pointlike approximation. Furthermore we 
expect an improvement from including a second atom- 
dimer scattering channel ~ Xad4>^ 4'''P''ip, cf. [l^. These 
improvements are expected to change the quantitative 
values of Od, ac2, Po, and sq, but not the qualitative sit- 
uation. We believe that already the present truncation 
will yield a reliable picture of the qualitative properties of 
the phase diagram, once it is extended to nonzero density 
and temperature. 
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